Polytypic nanowhiskers: electronic properties in the vicinity of the band edges 
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The increasing interest of nanowhiskers for technological applications has led to the observation 
of the zinc-blend/wurtzite polytypism. Polytypic nanowhiskers could also play, by their character- 
istics, an important role on the design of optical and electronic devices. In this work we propose 
a theoretical model to calculate the electronic properties of polytypic zinc-blend/wurtzite structure 
in the vicinity of the band edges. Our formulation is based on the k ■ p method connecting the 
irreducible representations in the interface of the two different crystalline phases by group theory 
arguments. Analyzing the composition of the states in the V point and the overlap integrals of the 
envelope functions we predict energy transitions that agree with experimental photoluminescence 
spectra. 



c 

o 

(N 
> 

(N 
(N 

o 



X 



INTRODUCTION 

Recently an increasing number of papers have re- 
ported high quality and control growth of nanowires and 
nanowhiskers (NWs) of III-V compounds such as InP, 
InAs, GaAs and GaP [H-0|. Nanoelectronics, nanopho- 
tonics and solar cells are some of the technological ap- 
plications for these nanostructures like described in the 
references above. 

Vapor-liquid-solid (VLS) method @ was used, in 1964, 
by Wagner and Ellis to grew a silicon whisker. Since then, 
other methods like molecular beam epitaxy (MBE) and 
metal organic vapor phase epitaxy (MOVPE) have also 
been applied to the growth of NWs. These techniques al- 
low the fabrication of homo and heterostructures in both 
axial and radial directions opening new opportunities to 
engineering novel device applications. 

A heterostructure is composed by assembling different 
materials side-by-side along a growth direction. The ho- 
mostructure is a similar concept, however, it is not the 
kind of material that is different but its crystalline struc- 
ture. This effect is also known as polytypism and have 
been observed in NWs structures of several III-V com- 
pounds. 

Under standard temperature and pressure conditions, 
the III-V compounds have as the bulk stable phase the 
zinc-blend (ZB) structure, exception made for the ni- 
trides that have the wurtzite (WZ) structure as the most 
stable phase. However, reducing the dimensions to the 
nanoscale level such as in these NWs, the WZ phase has 
been observed for non-nitride III-V compounds. The nat- 
ural explanation to this behavior is the small formation 
energy difference between the two crystal structures. 

The phase alternation creates ZB and WZ regions 
along the growth direction affecting the NWs optical and 
electronic properties, generating new forms of band off- 
sets and electronic minibands that can be explored to 
produce new kinds of applications. 

The WZ/ZB polytypism in III-V compounds can make 
a type II homostructure characterized by a staggered 
band-edge lineup, i. e., the charge carriers are spatially 



separated: the electrons will be concentrated in the ZB 
region while the holes will be in the WZ region. A sum- 
mary of NWs growth, properties and applications was 
made by Dubrovskii et al. [9j. 

From the theoretical point of view, the literature about 
these III-V NW structures, is focused in growth stabil- 
ity, surface energies and formation of polytypism (lp| - [l6| . 
Only recently theoretical calculations of the electronic 
properties in the band edges have been done for such 
structures (T^]. None has been done in the vicinity of 
the band edges. 

The polytypism was investigated in a paper of Mu- 
rayama and Nakayama [3] by ab initio calculations of 
the band offsets in WZ/ZB interfaces. A useful schemat- 
ics of the correspondence between the symmetry of the 
states of the WZ and ZB structures has been presented 
in the same article. 

The aim of this study is to propose a theoretical model 
based on group theory arguments and the k -p method to 
calculate the electronic band structure in the vicinity of 
the band edges and envelope functions of a WZ/ZB/ WZ 
homostructure for indium phosphide (InP). The effective 
mass parameters for the WZ bulk phase of InP were ex- 
tracted from a paper by Pryor and De that predicted 
the band structures of some III-V bulk compounds in the 
WZ phase (phosphides, arsenides and antimonides) . The 
ZB parameters are extracted from Vurgaftman et al. [2fj| . 

The formulation of a theoretical model for this new 
kind of nanostructures can lead to band gap engineering 
and novel electronic a and also work as a guide for the 
crystal growth community. 



THEORETICAL MODEL 

The existence two different crystal phases in a single 
whisker, can be explained by a simple joint analysis of the 
ZB and WZ crystal structures. Looking into ZB along the 
[111] direction and to WZ along the [0001] direction, one 
can found a noticeable similarity. Both can be described 
as stacked hexagonal layers. In fact, the ZB structure has 
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Figure 1. F point symmetry for the two crystal structures. 
Fhe numbers in parentheses are the degeneracy of the irre- 
ducible representation. Fhis figure was adapted form [l8| |. 



three layers in the stacking sequence (ABCABC) while 
the WZ structure has only two (ABAB). So, a stacking 
fault in WZ becomes a single segment of ZB and two 
sequential twin planes in ZB form a single segment of 
WZ 0]. 

An important difference between the two structures is 
the crystal field energy present in the WZ structure that 
defines the top of the valence band, neglecting spin, as 
a composition of two different bands, a non-degenerate 
and a bidegenerate bands. In ZB, the top of valence band 
is composed by a three-degenerate band. This difference 
in the energy bands between the two structures can be 
seen schematically for the T point in figure [TJ with and 
without the spin-orbit interaction. 

To calculate the band structure of a polytypic NW 
based in the k ■ p matrix formalism, it is necessary to un- 
derstand the differences between the two crystal struc- 
tures in the reciprocal lattice or, more precisely, what 
happens with the high-symmetry points in the ZB first 
Brillouin zone when we are dealing with the [111] growth 
direction and how they can be related to the WZ ones. 

Looking through the [111] axis in figure one can 
see that the L point of ZB structure is projected in the 
r point. Since the number of atoms in the primitive 
unit cell for WZ is twice the number in ZB, there are 
more states in the T point of WZ. This correspondence 
is done by assigning T and L points of the ZB to the 
r point of WZ. This compatibility can be seen in the 
energy diagram presented by Murayama and Nakayama 
18], that relates the irreducible representations of the ZB 
r and L points to the WZ T point. This energy diagram 
provides the essential insight for our method since we are 
going to relate k ■ p matrices developed with irreducible 
representations belonging to different crystal structures. 

The usual k-p matrix formulation is constructed in the 
vicinity of the T point for ZB and WZ structures 2l[ \22\ 
considering the first three valence bands and the first 
conduction band, represented by a 4 x 4 matrix. When 
the spin-orbit interaction is included this number is mul- 
tiplied by 2, resulting in a 8 x 8 matrix. These bands 
are exactly the ones presented in figure [TJ and represent 




Figure 2. (a) ZB and (b) WZ first Brillouin zones. The ar- 
rows (red online) represent the growth directions, (c) ZB 
conventional unit cell with two different coordinate systems, 
(d) WZ conventional unit cell with its common coordinate 
system. The [111] growth direction for ZB structure passes 
along the main diagonal of the cube and is represented in the 
primed coordinate system. 



a smaller set of the states first presented by Murayama 
and Nakamura 18]. Analyzing this, one may notice that, 
in the present method the irreducible representations of 
ZB [111] belong only to the T point and no irreducible 
representation of the L point will appear in the formula- 
tion. 

Another subject, asides the symmetry correspondence, 
must be pointed out: the k ■ p matrix for the ZB struc- 
ture is developed for the cubic main axes, which are the 
[100], [010] or [001] growth direction. This implies that 
the definition of the x, y, z coordinate system defines the 
matrix representations of the symmetry operations used 
to calculate the matrix elements in the k ■ p formulation. 
Changing the coordinate system to x' , y', z', which is the 
coordinate system of the WZ Hamiltonian, will enable us 
to relate the ZB and WZ matrices. 

We can either recalculate the matrix representations of 
the symmetry operations and develop the ZB k -p matrix 
or use the rotation matrices presented in the paper of 
Park and Chuang [23| to transform the ZB k ■ p Hamil- 
tonian in the x, y, z to the x' , y' , z' coordinate system. 

The basis, {|q)}, in the prime coordinate system is 
defined as (from this point on, the prime will be dropped 
out of the notation): 
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which is convenient to describe the WZ Hamiltonian. 

Since we are not treating the interband interaction ex- 
plicit in the matrix, the conduction band have a parabolic 
form for both crystals and is not affected by the rotation, 
reading as 



E c {k)=E g + E 
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where E g is the gap energy, Ey b m is the energy of the 
top valence band and for the ZB structure the electron 
effective masses are the same mjl = mjr- 

Applying the rotation matrices to the valence band 
ZB Hamiltonian we obtain a common matrix for the two 
structures: 
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The matrix terms arc defined as: 
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where fc-i 



v ± — n> x _i_ i,rvy. 

The only difference from the usual WZ formulation is 
the inclusion of a new parameter, A z , that takes into 
account the reorientation of the ZB growth axis. 
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Figure 3. (a) Bulk ZB valence band structure along the [111] 
direction. Jb) Bulk WZ valence band structure reproducing 
results in [191 ] along the [0001] growth direction. Both valence 
band structures were plotted for 20% of the first Brillouin 
Zone in k x and k z directions. Without spin-orbit the (x, y, z) 
coordinates belong to the same irreducible representation for 
ZB structure and then we can see that the band structure is 
isotropic in k x and k z directions. In WZ structure (x,y) and 
(z) coordinates belong to different irreducible representations 
making the band structure anisotropic in k x and fe z directions. 



The Ai(i = l,...,6,z), and A,(i = 1,2,3) WZ param- 
eters can be related to the — 1,2,3) and Ago ZB 
parameters comparing the ZB rotated matrix to the WZ 
matrix. 

In this work we will assume the WZ as ideal and there- 
fore the lattice parameters will considered the same as 
ZB. As a consequence, strain effects will not be taken 
into account. 



RESULTS AND DISCUSSION 

We first apply our model for bulk ZB and WZ InP to 
check the reliability of the matrix for both crystal struc- 
tures. This initial investigation is useful to fit the ef- 
fective mass parameters and to understand the valence 
band states in the band edge since now we are dealing 
with different basis states for the ZB matrix. 
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Energy Set of states 1 Set of states 2 

E 1 |ci> |c 4 ) 

E 2 a\c 2 )+/3\c 6 ) f3\c :i ) +a\c 5 ) 

E'j /3\c 2 ) -a\c 6 ) -a\c 3 )+/3\cs) 

Table I. The valence band energies and the respective eigen- 
states at the T point. 

The ZB parameters were obtained from Ref. [2(| and 
the WZ parameters were derived using the effective 
masses presented in Ref. flQ^ . Since in our model there 
is no explicit interband interaction only the valence band 
is presented in figure [3] for the bulk form of both crystal 
structures. 

The usual identification of the bands were used for ZB 
while in WZ we had to analyze the composition of the 
states in the T point. Table Q] summarizes the energies 
and states compositions in the T point that helps us to 
identify the WZ states. 



The energies and the a and (3 coefficients are given by 




Ei = Ai + A 2 (5) 



Energies a j3 

ZB E 1 =E 2 > E 3 0.57735 0.81650 
WZ Ei > E 2 > E 3 0.98345 0.18116 

Table II. The energy ordering and the a and f3 coefficients. 
The knowledge of these coefficients enable the identification 
of the states. 



ZB identification WZ identification 

HH HH 

LH CH 
SO LH 



Table III. The WZ identification for the ZB states. 



Since in our model we are using the WZ basis to de- 
scribe the ZB structure it is useful to think that ZB is 
a WZ structure with zero crystal-field energy. Noting 
this feature we can identify the ZB valence band states 
in the T point using the WZ identification. This helps 
us to understand the band edge in the WZ/ZB interface 
in a quantum well structure. Table IIIII relates the two 
different identifications. 
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The states identification were made as described in 
Ref. 24[: HH (heavy hole) states are composed only by 
\ci) or \ci), LH (light hole) states are composed mainly 
of \c2j or |cb) and CH (crystal-field split-off hole) are 
composed mainly of \cs) or \cq). According to the identi- 
fication of the energies presented in table HI E\ is the HH 
state and, if a > f), E2 is the LH state and E3 is the CH 
state or E2 is the CH state and E% is the LH state oth- 
erwise. The valence band ordering in the T point, which 
depends on the values of Ai, A 2 and A3, and the value of 
the coefficients for both crystal structures can be found 
in table HI 



As the dimensions of NWs are big enough to avoid lat- 
eral confinement, this structure can be considered as a 
2-dimensional system, like a quantum well or a superlat- 
tice. In this context, we simulated a single WZ/ZB/WZ 
polytypic quantum well with 160A well width and 320A 
barrier width so that the whole system has 480A. 

Using the band mismatch value given in Ref. [l9| it is 
possible to construct the quantum well profile. To un- 
derstand the confinement profile we plotted, in figure H^,, 
the band edge energy in the WZ/ZB interface. Since the 
composition of the states LH and CH is different in both 
crystal structures we did not think it was convenient to 
draw a line connecting the energies. We believe that 
there's a mix of these values to form the effective poten- 
tial. However, we can note that the conduction band and 
the first valence band form a type-II homostructure. 

It is also convenient to plot the diagonal terms of the 
total Hamiltonian to check the variation of the param- 
eters along the growth direction, which can be found in 
figure HJd. Although the nondiagonal terms of the matrix 
will mix the three profiles we can expect that the confine- 
ment of the electrons will be in the ZB region while the 
holes will be in the WZ region belonging in their most to 
the HH state. 

The band structure of the quantum well plotted up to 
10% of the k x direction and 100% of the k z direction is 
shown in figure [5] We can observe that the bound states 
in the valence band are mainly HH states, which means 
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Figure 4. (a) The band edge energy in the WZ/ZB interface 
relating the WZ and ZB states, (b) The diagonal potential 
profile of the total Hamiltonian along the growth direction. 



they are composed entirely by the |ci) and \a) states of 
the basis that have only \X) and \Y) components. This 
feature directly reflects in the luminescence properties of 
the system: the light will be almost fully polarized in 
the xy plane, i. e., transversal to the growth axis of the 
polytypic quantum well. The LH states have less than 
5% of \Z) component. 

The probability densities calculated from the envelope 
functions for the L point of the ELI, EL2, EL3, HH1, 
HH2 and HH3 bands are shown in figure H>] These results 
agrees with the experimental data suggesting a spatial 
separation of the carriers, due to the fact that is more 
likely to find holes in the WZ region and electrons in the 
ZB region. 

In a quantum well the selection rules for dipole optical 
transition are given by the matrix element M 25 1. 
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M = (E\r- e\ H) = M c M g (10) 

This matrix element can be divided in two parts: M c 
and M g . M c is the matrix element between the electron 
and hole states in the basis {|q)} that depends on the 
polarization of the light and M g is the overlap integral 
between the electron and hole envelope functions. 



\r ■ e c h ) 
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M g {%) Energy difference (meV) 



EL1-HH1 
EL1-HH3 
EL2-HH2 
EL3-HH1 
EL3-HH3 



-2.6 
-6.4 
~12 
~23 
-33.5 



1370.76 
1372.89 
1405.06 
1453.86 
1455.99 



Table IV. Finite overlap integrals (M g ) and energy difference 
between some conduction and valence band states at the T 
point. We can note that overlap is bigger for excited states. 



M g = / g* e (z)g h (z)dz (12) 

The M c matrix elements are known from symmetry 
properties of the L point and the non-vanishing elements 
are (S\x\X), (S\y\Y) and (S\z\Z). We just have to 
calculate the overlap integral M g . Since the potential 
profile is even, the envelope functions have well defined 



parities, they are either even or odd and we can expect 
to have null overlap integrals between conduction and 
valence states. The non-zero overlap integrals and the 
difference energy in the L point are presented in table 
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Figure 6. The probability densities of the electrons (a) and 
holes (b). We can note the spatial separation of the carriers. 



Although we have achieved reliable results with the 
model considered here for InP, we believe that improve- 
ments should be made to treat other III-V compounds. 
In systems where the lower conduction band has symme- 
try Ti in the WZ phase and the spin-orbit splitting is 
bigger or equivalent to the gap energy (InAs and InSb 
fljj p the full 8x8 k-p method should be considered with 
the inclusion of the Kane interband momentum matrix 
elements (P) and also Kane quadratic terms (B). In sys- 
tems with symmetry in the lower conduction band 
for WZ phase, like AlAs or GaAs [l9|, the full symmetry 
analysis should be redone. 



CONCLUSIONS 

In summary, we have presented a theoretical method 
based in group theory arguments and the k ■ p method 
to calculate the electronic properties of InP polytypic 
structures in the vicinity of the band edges. 

Understanding the behavior of the band structure of 
the calculated system we could predict the trends in the 
luminescence energy transitions. The preliminary results 
are in good agreement with the experimental data avail- 
able so far such as carriers' spatial separation and light 
polarization in luminescence spectra. 

With the proposed model, it will be possible to predict 
the optical properties of selected structures and to design 
systems that best fit new device applications. 
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